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PREFACE 


Slender  re-entry  vehicles  at  high  altitude  and  high  velocity  pos¬ 
sess  a  flow  field  which  results  from  the  interaction  between  the  boundary 
layer  formed  near  the  surface  and  the  shock  wave  originating  from  the 
leading  edge  or  tip.  This  Memorandum  studies  the  analytical  conditions 
necessary  for  the  solution  of  a  set  of  simplified  equations  which  describe 
the  interaction  flow  field.  Surface  mass  transfer,  which  occurs  in  the 
case  of  an  ablating  vehicle,  is  also  considered.  The  results  of  this 
Memorandum  wilT.  be  useful  in  the  interpretation  of  experimental  data 
and  the  implementation  of  numerical  analyses.  The  Memorandum  is  part 
of  a  continuing  study  for  the  Advanced  Research  Projects  Agency  in  re¬ 
entry  aerodynamics . 

One  of  the  authors,  T.  Y.  Li,  is  a  professor  in  the  Department  of 
Aeronautical  and  Astrorautlcal  Engineering,  Ohio  State  University, 
Columbus,  Ohio,  and  consultant  to  The  RAllD  Corporation. 
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ABSTRACT 


An  analytical  formulation  of  the  problem  of  hypersonic  viscous 
interaction  on  a  very  slender  body  of  revolution  with  a  thick  boundary 
layer  and  mass  transfer  is  attempted.  Particular  attention  is  directed 
to  the  establishment  of  the  mathematical  restrictions  necessary  to  ensure 
the  existence  of  similar  solutions  of  the  laminar  boundary  layer  equa¬ 
tions  for  this  class  of  problems.  Yasuhara's  analysis  is  extended  to 
include  surface  mass  transfer  and  binary  mixture  effects,  and  it  is 
shown  that  similar  solutions  are  possible  only  for  three-quarter-power- 
law  bodies  of  revolution.  With  regard  to  the  problem  of  hypersonic 
vis  ous  interaction  on  a  slender  body  of  revolution  with  a  thick  boundary 
layer  and  surface  mass  transfer,  it  is  shown  that  the  criterion  for 
the  existence  of  siiuilar  solutions  is  that  t^v^  x^  must  govern  the  body 
surface  shape  and  the  surface  injection  velocity. 
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SYMBOLS 

“  constant  defined  in  Eq.  (72) 

=  constant  defined  in  Eq.  (77) 

=  constant  defined  in  Eqs.  (86)  and  (133) 

»  constant  defined  in  Eqs.  (118)  and  (134) 

■  constant  defined  in  Eq.  (81) 

'  dimensional  constant  defined  in  Eq.  (96) 
concentration 

specific  heat  at  constant  pressure 

constant  defined  in  Eq.  (119) 

coastant  of  integration  defined  in  Eq.  (127) 

diffusion  coefficient 

constant  defined  in  Eq.  (120) 
dimensionless  enthalpy  >=  H/He 
total  enthalpy 
specific  enthalpy 

constant  defined  by  definite  integral  given  in 
effective  number  of  degrees  of  freedom 
dimensionless  velocity  in  x-direction  =  u/u 
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thermal  conductivity 
characteristic  dimension 
Lewis  number 
Mach  number 

constant  exponent  defined  in  Eq.  (96) 
molecular  weight 

Prandtl  number 
pressure 

universal  gas  constant 
mixture  gas  constant 
radial  distance 
Schmidt  number 
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similarity  coordinate  defined  in  Eq.  (12) 
similarity  coordinate  defined  in  Eq.  (10) 

similarity  coordinate,  defined  in  Eq.  (14) 
flow  velocity  component  in  x-direction 
flow  velocity  component  in  y-direction 

coordinates  as  shown  in  Fig.  1 

dimensionless  concentration  C. /C. 

1  iw 

edge  of  the  boundary  layer 
i  component 
generalised  coordinate 
wall 

undisturbed  stream 
angle  defined  in  Fig.  1 

pressure  gradient  parameter  defined  in  Eq.  (33) 
specific  heat  ratio  C  /C 

p  V 

boundary  layer  thickness 

displacement  thickness  of  boundary  layer 
hypersonic  parameter  defined  in  Eq.  (71) 
dimensionless  similarity  coordinate  defined  in  Eq.  (13) 
dimensionless  similarity  coordinate  defined  in  Eq.  (11) 

dimensional  similarity  coordinate  defined  in  Eq.  (15) 
constant  defined  in  Eq.  (72) 

function  defined  in  Eq.  (24) 

viscosity 

density 

stream  function 

parameter  defined  in  Eq.  (143) 
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I.  INTRODXTION 


The  purpose  of  this  Memorandum  is  to  present  an  analytical  formu¬ 
lation  of  the  problem  of  hypersonic  viscous  interaction  on  a  very  slen¬ 
der  body  of  revolution  with  a  thick  boundary  layer  and  surface  mass 
transfer.  Particular  attention  is  directed  to  establishing  the  mathe¬ 
matical  restrictions  necessary  to  ensure  the  existence  of  similar  so¬ 
lutions  of  the  laminar  boundary  layer  equation',  for  this  class  of  prob¬ 
lems.  The  determination  of  similarity  conditions  is  important  because 
similar  solutions  are  amenable  to  Cv.aparatively  simple  anaiysis^^^  and 

find  diverse  application  through  the  use  of  the  local  similarity  con- 

(2) 

cept.  Recently,  Yasuhara  demonstrated  chat  similar  solutions  can 
be  determined  for  a  slender,  three-quarter-power- law  body  of  revolu¬ 
tion  when  the  boundary  layer  thickness  is  of  the  same  order  as  the  body 
radius.  In  this  Memorandum  we  extend  Yasuhara 's  analysis  to  include 
surface  mass  transfer  and  binary  mixture  effects,  and  show  that  in  the 
extended  case  similar  solutions  are  possible  only  for  three-quarter- 
power-lav  bodies  of  revolution.  It  is  further  shown  that  these  simi¬ 
lar  solutions  are  characterized  by  the  restrictive  condition  that 
,  * 

r^/r^  =  constant.  This  requirement  for  similar  solutions  is  particu¬ 
larly  pertinent  in  the  regime  where  r^/r^  ■  0(1)  or  '  0(1). 

The  equations  which  describe  the  boundary- layer  growth  on  a  body 

2  2 

of  revolution  contain  the  transverse  curvature  expression  (r  /r^)* 

For  similar  solutions  to  exist,  this  expression  must  be  a  function  of 
the  similarity  variable  T).  Particular  attention  is  directed  to  the 
regime  where  «  1  or  Yg/r^  !•  We  derive  a  new  system  of  equa¬ 

tions  that  possesses  the  same  analytical  features  as  the  equations 
studied  previously  by  Stewartson, Glauert-Lighthill, and  Mark.^^^ 
The  possibility  that  this  new  system  has  similar  solutions  is  then  ex¬ 
amined  under  appropriate  boundary  conditions,  including  the  injection 
of  gas  at  the  body  surface.  The  results  indicate  that  similar  solu¬ 
tions  characterized  by  prv-similarity  in  the  boundary  layer  can  be  ob¬ 
tained.  Finally,  with  regard  to  the  problem  of  hypersoni.c  interaction 


See  List  of  Symbols 
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on  a  very  slender  body  of  revolutl.on  with  a  thick  boundary  layer  and 

surface  mass  transfer,  it  is  shown  that  the  criterion  for  the  exist* 

1/2 

ence  of  similar  solutions  is  that  ~  x  '  must  govern  the  body  sur¬ 

face  sh.ipe  and  the  f.  irface  injection  velocity. 
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II.  CLASSIFICATION  OF  FLOW  REGIMES 


Viscous  flow  regimes  may  be  usefully  classified  with  respect  to 
the  ratio  of  the  characteristic  body  radius  to  the  thickness  of  the 
boundary  layer.  Three  such  regimes  described  below,  may  be  distin¬ 
guished  for  the  purposes  of  our  problem. 


THE  M  REGIME 

This  regime  cotaprises  cases  in  which  ”  '^(e)  where  e  «  1. 

Here,  Mangier '&  transformation^^^  permits  a  reduction  of  axi- 

symmetrij  flow  boundary  layer  equations  to  two-dimensional  boundary 
layer  equations.  Transverse  curvature  effects  are,  of  course,  unim¬ 
portant  in  this  regime. 

THE  P-E-Y  REGIME 

This  regime  comprises  cases  in  which  y  /r  “  0(1)  or  r  /r  ■  0(1). 

e  w  e  w 

The  appropriate  compressible  boundary  layer  equations  with  mass  trans¬ 
fer  are  considered  by  Yasuhara^^^  and  Probstein  and  Elliott, who 
show  the  importance  of  the  effect  of  transverse  curvature  in  this  re¬ 
gime. 

THE  G-L  REGIME 

In  this  regime  che  boundary  layer  thickness  is  large  compared  with 

the  radius  of  the  body  of  revolution;  that  is,  y  /r  »  1.  Using  the 

e  w  _ 

criterion  stated  above,  we  can  characterize  this  regime  by  = 

0(e),  where  e  «  1.  Stewartson^^^  and  Glauert  an'’  Lighthill^^^deait 
with  incompressible  flow  problems  in  this  regime.  Mark^^^  and  Steiger 
and  Bloom^  also  considered  certain  cic-ses  of  compressible  boundary 
layer  flows  in  this  regime.  Their  results  demonstrated  the  extreme 
importance  of  transverse  curvature  effects.  The  boundary  layer  equa¬ 
tions  for  the  G-L  regime  possess  a  basic  structure  which  leads  to  a 
logarithmic  velocity  profile  distribution  near  the  body  surface  in 
the  case  of  zero  mass  transfer. 
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III.  BINARY  BOUNDARY  LAYER  EQUATIONS  AND  BOUNDARY  CONDITIONS 

The  no-slip  condition  must  apply  at  the  body  surface..  This  re¬ 
quires  the  use  of  the  boundary  layer  equations  in  the  rtuion  adjacent 
to  the  body  surface.  A  gas  is  injected  normal  to  the  surface  such  that 
a  binary  boundary  layer  flow  results.  The  binary  boundary  layer  equa¬ 
tions  for  axisjinnetric  flow  may  be  written: 


where 


H 

b 


specific  heat  at  constant  pressure  of  the  i  component 
E  C.C 

1  p. 

1 

mass  fraction  of  i^^  component 
binary  diffusion  coefficient 

2 

total  enthalpy  =  h  +  (l/2)u 

Z  C,h. 

1  1 
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til 

®  enthalpy  the  i*"  component 

k  =  thermal  conductivity  of  the  fluid 
Le  =  Lewis  number*  p  C  /k 

12  p 

Fr  -»  Prandcl  number,  n  C^/k 

p  »  fluid  Pressure 

r  =  r  +  V  crs  a 

w  ^ 

defines  the  shape  of  the  body  surface 

Sc  =  Schmidt  number,  M./p  Dj^2 

u,v  *  flow  velocities  in  the  x,y  directions,  respectively 
x,y  =  coordinates,  as  shown  in  Fig.  1 
p  =  viscosity  coefficient  of  th';  fluid 
p  =  fluid  density 


Within  the  boundary  layer  we  have,  from  Eq.  (3), 


Pg(x)  »  RpT 


(6) 


where 

(7) 

m^  =  molecular  weight  of  i-gas 
R  =  universal  gas  constant 

R^  =  component  gas  constant 
R  =  mixture  gas  constant 

The  boundary  conditions  at  the  wall  are  characterized  by  the  no¬ 
slip  condition,  constant  wall  temperature,  the  concentration  of  the  in 
jectant  gas,  and  the  Eckert  condition,  which  specif’ es  that  the  normal 
mass  flow  velocity  of  the  main  component  (i  =  2)  in  the  gas  mixture 
vanishes  at  the  body  surface. 


R  »  E  C.R. 

1  1 

1 


E  C  — 
^  i  ra. 


-6- 


when  y  ■  0 ;  u  ■  0 


C,  ■  C,  “  constant 
1  Iw 


(8a) 

(8b) 


h  ■  H  ■  constant 
w 


(8c) 


(pv) 


w 


C  '*12  W  J 

zw 


(8d) 


w 


At  the  edge  o£  the  boundary  layer,  the  properties  of  the  flow  must 
match  those  given  in  the  inviscid  flow. 


when 


6. 


(9a) 


"l  =  ° 


H  -  H 


(9b) 

(9c) 


TRANSFORMATION  OF  BINARY  BOUNDARY  LAYER  EQUATIONS 

We  shall  introduce,  a  coordinate  transformation  to  translorm,  under 
certain  prescribed  conditions,  the  partial  differential  equations  de¬ 
scribing  the  axially  symmetric  laminar  boundary  layer  into  a  set  of 
nonlinear  coupled  ordinary  differential  equations.  The  transformation 
is  quite  general  so  that  it  can  be  applied  to  all  the  regimes  described 
above.  The  boundary  layer  coordinates  (x,y)  are  transformed  into  the 
similarity  coordinates  follows: 


=  r  Cp  u  p.  rf  dx 
k  Jq  e  e  e  k 


r  f-  r  dy 


k  Jq  pg 


(10) 


(11) 
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where  C  =  the  Chapraan-Rubetin  constant 
“  characteristic  dimension 

The  transformation  can  be  applied  to  the  P-E-Y  and  the  G-L  regimes  as 
follows : 

In  the  P-E-Y  regime,  r.  -  r  ,  thus 

K  W 


Cp  p  u  r  dx 
e  e  e  w 


P  u  ^ 


In  the  G-L  regime,  r^^  «  r^,  thus 


(12) 


(13) 


f'  Cp  M.  u  r^ 


(14) 


P  u 
e  e 


P 


«  *  Jn  P 
2s  0  '^e 


dy 


(15) 


We  Shall  discuas  late,  these  fonss  o£  the  tra„s£et™tlo„  and  their  slg. 
nificance.  Their  immediate  application 


tions  which  can  now  be  rewritten  in  the 
stream  function  t  is  defined  by; 


is  to  the  boundary  layer  equa- 
new  variables  (s^,Tlj^).  The 


dx  = 


(16) 


S'!; 

dy 


Pur 


(17) 


We  consider  those  flows  where  n  h  r  , 

s  wnere  u,  H,  and  are  expressed  £s  follows! 


t  '=k> 


(18) 


^0. 


■J 


8- 


H-  “ 

e 


cZ  = 


where  the  prime  denotes  Equations  (2),  (4),  and  (5)  may  now  be 

written: 


f  -r^  \  2s,  du  ,  -  p  . 

(x  ~  K")  +  KK"  =  — — £  (k'2  .  _e^ 

V  2  /  V  p  , 


e  r„  ' 


(x^|-Z^')  +K2'  .  2d^(|i  K'  -  Zj') 


where 


-  i  PM- 

C  p^M- 
e  e 


\ 


(24) 
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In  the  new  coordinate  system,  the  boundary  conditions 
follows : 


are  expressed  as 


when  "H  =  Q 
k 


K'(0)  =  0 

(25) 

2^(0)  -  1 

(26a) 

H 

G(0)  = 

^  ^  H 

e 

(26b) 

2  s. 


dK 


when  =  n 
k  'k 


(27) 

K"(T1  )  =  1 

e 

(28) 

)  -  0 
e 

(29a) 

G(\  )  »  1 

e 

(29b) 

10. 


The  lowing  relationships  were  used  in  the  derivation  of  Eqs.  (34) 
and  (35) : 


11- 


R 

R 


1  + 


^lw^l^\’®k^ 


(36) 


ji  +  2 


(37) 


The  j.  is  defined  as  the  effective  number  of  degrees  of  freedom  of  the 
1  component.^  In  this  study,  j  =  5  for  diatomic  gases  and  j  -  3 
for  monatomic  gases.  Substitution  of  Eqs.  (30)  to  (35)  into  Eqs.  (21) 
to  (23)  gives; 


2  s' 

^  K")  +  kK 


Ye  -  1 

2Y  ®jj) (G  -  ) 


(38) 


+  2s Tk'  |2-  -  |!L  G'1 

kL  as^  J 

^  &  -  0  ^2^\« 
k 


9  * 

—  —  z') 

2  Sc  ^iJ 


+  KZ, 


fdZi 

‘‘VST*'' 


(39) 

(40) 


For  the  problem  of  hypers 
revolution  with  injection 


onic.  viscous  flow  past  a  very  slender  body  of 
of  coolant  gas  at  the  surface ,  these  equations 
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can  be  treatea  under  the  boundary  conditions  given  in  Eqs.  (25)  to 
(29). 

If  similar  solutions  are  to  be  obtained,  the  following  conditions 
must  be  satisfied: 

(a)  K  - 

G  -  Gy^\) 


Hence,  Eqs.  (38)  to  (40)  become: 


+  KK" 


PFl(\)(G  -  K'2) 


(38a) 


k 


-  iV] 

(40a) 


These  equations  are  to  be  solved  under  the  boundary  conditions  given 
in  Eqs.  (25)  to  (29),  whence.  Eq.  (27)  becomes 


K(0) 


(27a) 
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(b) 

X  »  X(\) 

..1) 

Sc  -  Sc(71j^) 

(42) 

Pr  -  Pr(Tlj^) 

(43) 

(c)  The  transverse  curvature  parameter  r^/r^  must  be  a  function  of 
Tlj^;  that  is, 

h  *  (44) 

(d)  The  pressure  gradient  parameter,  as  defined  in  Eq.  (33),  must  be 
constant: 


d  Xn  p 

^  ^  d  Xn  s^  (45) 

k 

(e)  The  surface  boundary  conditions  are  defined  such  that  C  H  /H 
/  /  2  Iw*  e* 

snd  (r^/rj^)  are  constants. 

THE  ANALYTIC  CHARACTERISHCS  OF  THE  G-L  EOUArmNS 

Glauert  and  Lighthill^^^  studied  the  problem  of  the  incompressible 
axisyrametric  boundary  layer  of  a  long  cylinder.  Their  primary  purpose 
was  to  examine  a  very  thick  boundary  layer  on  a  slender  body;  that  is, 

where  rjr^  »  1.  in  our  notation,  their  equation  for  the  boundary 
layer  flow  is 


[TlV]'  +  iKK"  =  0 

Suppression  of  the  pressure  gradient  terra  in  Eq.  (38a)  and  taking 
X  =  1  gives 


(46) 


"14“ 


-  0 


C47) 


For  incompressible  flow,  it  will  be  shown  later  that  (r/r  =  T]*/!]*. 
Hence,  Eq.  (47)  can  be  written  ^  ® 


0 


(48) 


The  similarity  between  Lqs.  (46)  and  (48)  demonstrates  that  thi  system 

of  equations  presented  here  is  entirely  consistent  with  the  Giauert 
and  Lighthill  equation* 


SI.MILARITY  CONDITIONS  FOR  THE  NORMf*L  VELOCIiY 

!n  seeking  similar  selutlerm  Cor  problems  relating  to  mass  trans- 
£er  in  the  boundary  layer.  It  Is  essential  to  ensure  that  the  normal 
velocity  at  the  surface,  which  governs  the  axial  rate  of  mass  transfer 
Is  a  function  of  the  similarity  variable  only.  This  can  be  shown  by 
first  combining  Eqs.  (11)  and  (17); 


u 


(49) 


Combining  Eqs.  (18)  and  (49)  gives. 


in  the  case  of  similar  solutions. 


i'  =  K(71j^) 


(50) 


From  Eq.  (16),  one  obtains 


-pvr 


(51) 
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No« 


9y  as,  ax 

k 


(’2) 


Inversion  of  Eq.  (11)  gives 


y  = 


—ii 

Wk  Jo 


dTl 


(53) 


which  by  differentiation,  yields 


^  '  ^“e'k  ■  P  “  r.  sT  ^  ■i’l  (54) 


e  e  k  k 


Combining  Eqs.  (52)  and  (54)  yields 


5 

9x 


_ k  /,  _  ^  k  a  I  Dr  A  P 

2s.  ax  t  ”  p  u  r  Qs  7^ —  f  ^ 


dT|  (55) 


This  result  is  substituted  into  Eq.  (51)  to  yield 


Dvr  rr  ^®ir  ^(P_“  r,  ) 


/2  s,  Bx 
k 


■e-e^'k 


From  Eqs.  (30)  and  (33),  we  obtain,  respectively 
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Heuce,  Eq.  (56)  can  be  written 
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where 


(60) 
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In  hypersonic  boundary  layer-*  the  terms  of  0(1/M^)  are  neglected,  and 
EqSi  (59)  and  (60)  reduce  to 
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pur, 
e  e  k 


u  I,  ; 
e  e  k 


where 
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It  should  be  pointed  out  that  the  right-band  side  of  Eq.  (61)  is 
a  function  of  only,  provided  chat 
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Now  condition  (b)  can  be  written  as  follows; 
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as  defined  in  Eq.  (57).  In  the  hypersonic  adiabatic  free  stream,  the 
total  enclialpy  is  constant  so  that  dh^  +  du^  =  0.  Indeed,  in  the 
hypersonic  boundary  layer,  6  «  6*  so  that  very  little  mass  flov-s  into 
the  boundary  layer.  Hence,  the  specific  entropy  is  constant  along  the 
edge  of  the  boundary  layer  and  it  follows  that 
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(66) 


./hich  leads  to 
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We  therefore  have 
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This  shows  that 
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const. 


(68) 


(69) 


(70) 


is  sufficient  for  condition  (b) ;  this  is  a  geometrical  constraint  on 
the  problem. 

BOUNDARY  LAYER  EFFECTS  ON  THE  MATCHING  CONDITIONS  BETWEEN  INVISCID 
VISCOUS  FLOW  REGIONS 

« 

The  axis  of  symmetry  of  the  body  lies  along  the  direction  of  the 
undisturbed  free  stream.  At  the  body  surface,  a  parcicular  rate  of 
mass  transfer  due  to  -oolant  gas  injection  is  specified.  The  disturbed 
flow  region  surrounding  the  body  is  separated  from  the  undisturbed  flow 
by  the  leading-edge  shock  wave.  The  region  between  the  shock  wave  and 
the  body  surface  can  be  divided  into  an  inviscid  flow  region  and  a  vis¬ 
cous  boundary  layer  region.  We  shall  match  the  normal  velocity  and  the 
pressure  on  the  common  boundary  of  these  inviscid  and  viscous  flow  re¬ 
gions. 

Stewartson^^^^  and  Oguchi^^^^  have  obtained  similar  solutions  of 
the  hypersonic  inviscid  flow  region  using  the  small-perturbation  theory. 
These  results  can  be  summarized  in  the  following  equations  for  the  nor¬ 
mal  velocity  and  pressure  in  the  flow  field: 
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(71) 
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e  00 
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where  A,  N,  and  9  are  constants.  These  values  v  and  p  will  be 
o  e  e 

matched  with  the  corresponding  viscous  flow  solution  at  the  edge  of 
the  boundary  layer.  These  matching  conditions  will  determine  n  and  g 
in  Eqs.  (71)  and  (72). 

The  value  of  can  be  obtained  from  the  boundary  layer  analy¬ 

sis  by  evaluating  Eq.  (56)  at  y  =  6.  This  result  has  also  been  obtained 
by  the  authors  usie^j  a  direct  integration  of  the  continuity  equation, 
and  gives,  for  h3rpersonic  boundary  layers. 
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v,  SIMEIAR  SOLUTIONS  IN  THE  P-E-Y  REGIME 

We  note  that  iu  the  P-E-Y  regime,  =>  and  become 

(s,T)),  as  defined  in  Eqs«  (12)  and  (J’).  The  conditions  for  similar 
solutions  in  this  regime  are: 


X  -  m) 


(74a) 


Sc  =>  sccn) 


(74b) 


Pr  =>  PrCH) 


(74c) 


2  2 

(b)  The  transverse  curvature  parameter  r  /r  must  be  a  function  of  Tl: 

w 


(c)  The  pressure  gradient  parameter  as  defined  in  Eq.  (33)  must  be 


constant: 


d  ;6n  p 

_ ^ 

d  Xn  s 


(d)  The  surface  boundary  conditions  are  defined  such  that  and 

H  /H  are  constant.  This  requires  uniform  distribution  of  the 
we 

wall  temperature  and  the  concentration  of  the  injectant  gas. 

(e)  The  boundary  layer  effects  given  in  Eq.  (73)  must  be  included  in 
the  matching  condition. 

The  above  condition  (c)  implies  that 


P/2  * 

s  =  A  p 


where  A  is  a  dimensional  constant.  From  Eq.  (12),  one  obtains 


21- 


s 


r 


li_  'Is. 

RT  X 


2  ^ 

p  r  d.v 
e  w 


(78) 


In  the  hypersonic  flow  regime,  u  u  .  In  Eq.  (78),  X  can  be  treated 

Q  00 

as  a  constant  evaluated  at  some  reference  temperature  and  composition, 
X  =  X^.  To  the  same  approxiv  ation,  then. 


RT 


and  hence 
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where 
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const. 
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Combining  Eqs.  (77)  and  (80)  yields 

*  2/8 
(A  ) 

t -  =  const.  (82) 


Consider  now  condition  (b) ,  which  requires  that  r^/r^  be  a  func- 

w 

tion  of  T|  only.  From  geometry 
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1  .  (2  cos  Of) 
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(84) 


Equations  (58),  (80),  and  (84)  together  give 


(^)  “  1  +  (v^v/2  u^  cos  c/)  (-~^)  a*  F^(G  -  K'^)  dTl  (85) 


2  ^ 
where  terms  of  0(1/M^)  have  been  neglected,  cos  O'  =“  1,  and  a  is  de¬ 
fined  as 
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(86) 


*  2,2 

Now  a  must  be  a  constant  in  order  that  r  /r^  be  a  function  of  T)  only. 

Consider  finally  the  matching  condition,  from  Eqs.  (71)  and  (73), 

we  liave 
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In  matching  the  viscous  and  inviscid  solutions,  we  have  included  the 
boundary  layer  effects.  Since 


r  =  r  +  6  cos  O'  (88) 
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Equation  (87)  may  be  written 
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Since,  by  Eq.  (51), 
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we  may  write 
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where  B  =  (p,u/RTX)  .  The  final  expression  for  tiie  matching  condition 
is  therefore 


ne0 
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Similar  solutions  in  the  P  E-Y  regime  must  satisfy  the  conditions 
given  in  Eqs.  (82),  (86),  and  (93).  Equation  (72)  can  be  rewritten 
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Let 


r  =  bx 
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tn 


(96) 


where  b  is  a  dimensional  constant.  From  Eqs.  (82),  (86),  and  (93), 
one  obtains 
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2 
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(99) 

n  =  m  =  3/4,  8  =  -1/2. 

Equations  (38a),  (39a), 

and  (40a)  may  now  be  written 


(x  4 


Y  -  1 


+  KK  = 


4Y 


F^(K'^  -  G) 


(100) 


w 


+KC'  .  2  [x  4^7-  O'*"] 


(101) 


W 


w 


+>“1'  -  0 


w 


(102) 


The  boundary  conditions  for  this  set  are: 
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( 103a) 
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(104c) 


The  geometrical  constraint  for  similarity  given  by  Eq.  (70)  will 
now  be  examined  for  this  regime.  Equation  (80)  reduces  to 
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By  Eqs.  (94)  and  (96),  we  obtain 


dr 


w 


dx 


ds 

dx 


dr 

_ V 
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2(n  -  1)  +  2m  +  1 


=  const.  (106) 


Thus  the  condition  of  similarity  is  satisfied.  The  foregoing  arguments 
show  that  similar  solutions  exist  for  the  case  of  a  slender  body  of 
revolution  in  hypersonic  flow  with  surface  mass  transfer.  The  shape 
of  the  body  must  be 
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r  =  (107) 

w 

For  zero  mass  transfer,  the  result  confirms  the  conclusions  reached  by 
Yauuharar 
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VI.  SIMILAR  SOLUTIONS  IN  THE  G-L  REGIME 


In  the  G-L  regime,  r,  »  r  and  (s,  ,1]  )  io  transformed  into 
*  *  R.  e  k  k 

(s  ,Ti  ),  as  shown  in  Eqs.  (14)  and  (15).  The  conditions  for  similar- 

ity  in  this  regime  are: 

(a)  X  -  X(ll*) 

(108a) 

Sc  »  Sc  (11*) 

(108b) 

Pr  =  Pr(T1*) 

(108c) 

2  2 

(b)  The  transverse  curvature  parameter  r  /r 

e 

must  be  a  function  of  Tj 

(109) 
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(c)  The  pressure  gradient  parameter  as  defined  in  Eq.  (33)  must  be 
constant: 
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M)  The  viscous  term  given  in  Eq.  (73)  must  be  included  in  the  match¬ 
ing  condition: 
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(e)  The  surface  boundary  conditions  are  defined  such  that  C,  and 

Iw 

H  /H  are  constant, 
w  e 

(f)  The  boundary  condition  [Eq.  (27a)]  relating  to  the  mass  transfer 
at  the  surface  must  be  compatible  with  the  requirements  for  a  sim¬ 
ilar  solution. 
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Condition  (c)  implies  t’lat 
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where  Aj^  is  a  dimensional  constant.  Equations  (14),  (81),  and  (112) 
give 


dx 


.  2/8 

(Aj^y 


const. 


(113) 


which  is  similar  to  Fq.  (82)  for  the  P-E-Y  regime. 

Condition  (b)  requires  that  the  transverse  curvature  parameter 
2  2  ^ 

r  /r^  must  be  a  function  of  T]  .  Since  this  study  concerns  a  slender 
body  of  revolution,  an  approximation  is  made  by  assuming  that  the  x- 
axis  and  the  axit  of  symmetry  are  coincident;  that  is,  cos  o  =“  1. 
Equation  (15)  then  yields 
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where  denotes  the  boundary  layer  thickness  in  the  T]  variable. 
These  equations  show  that,  in  the  G-L  regime. 


Combining  Eqs.  (58)  and  (115)  gives 
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2 

where  terms  of  0(1/M  )  are  neglected.  Equation  (116)  shows  that 
2  ®  * 

(••'r^)  is  a  function  of  T]  within  the  present  approximation.  Indei  d, 
r/r^  —  y/6,  so  that  r/r^  is  a  proper  similari"y  variable  in  the  physi¬ 
cal  coordinate  of  the  boundary  layer. 

For  this  regime,  the  matching  condition  in  Eq,  (87)  is 
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Thereforei  similar  solutions  in  the  G-L  regime  must  satisfy  the 
conditions  in  Eqs.  (113),  (117),  and  (118)«  Next  we  assume  that 
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cx 


(119) 
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e 


dx 
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and  determine  that  m  ■  n  *  3/4,  3  »  -  1/2.  The  boundary  layer  equa¬ 
tions  may  be  rewritten  as: 
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The  boundary  conditions  are 

=  0  K'(0)  =  0  (124a) 

Zl(0)  =  1  (124b) 
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and 
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Condition  (f)  requires  that 
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const. 
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(125b) 
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Treatment  of  the  geometrical  constraint  given  by  (70)  gives 
the  same  result  for  the  G-L  regime  that  was  shown  previously  in  Eq. 
(106).  In  the  G-L  regime,  therefore,  '‘exact''  similar  solutions  are 
obtainable  only  for  the  case  where  r^  ~  For  very  slender  bodies, 

the  condition  (r^/r^)  «  0  should  be  considered.  This  suggests  an  in¬ 

consistency,  Since  K(0)  =  0  for  finite  Z/(0).  To  alleviate  this  diffi- 
culty  the  present  authors,  in  Ref.  18,  adopt  the  following  condition: 


K(0) 


(127) 


where  D  is  a  finite  and  negative  constant.  When  this  alternative  con 
dition  can  be  accepted,  it  replaces  Eq.  (124d)  thus  raakiivg  then  pos¬ 
sible  "approximate"  similar  solutions  of  the  G-L  regime.  Within  this 
approximation,  it  may  be  observed  that 
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in  the  hypersonic  viscous  interaction  solutions  of  the  G-L  regime. 

This  approximate  result  is  completely  analogous  to  the  strong  inter¬ 
action  solution  on  a  hypersonic  flat  with  surface  mass  trans- 

(14) 

ter.  Pressure  data  on  a  solid  cone  in  the  st: ong  interaction  re¬ 

gion  have  been  analyzed  by  Yasuhara, and  these  data  are  in  general 
agreement  with  the  present  prediction.  This  agreement  may  be  fortui¬ 
tous,  however,  as  the  hypersonic  viscous  interaction  problem  for  a 
slender  impermeable  cone  yields  rigorously  only  a  nonsimilar  solu¬ 
tion. 
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VII.  MATCHING  THE  INVISCID  AND  VISCOUS  FLOW  SOLUTIONS 


Consideration  of  the  above  solutions  and  those  of  Mirels^^^^  for 
viscous  and  inviscid  regions  show  that  similar  solutions  exist  both 
m  the  inviscid  and  viscous  flow  regions.  It  remains  now  to  match 
uhese  solutions  to  obtain  the  proper  behavior  in  the  strong  interac¬ 
tion  region. 

Equations  (83)  and  (85)  yield 
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In  the  P-E-Y  regime,  \  =  I],  a,  =  a*,  r,  =  r  . 
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In  the  G-L  regime, 
shown  that 
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b,  c,  and  d  being  the  same  as  defined  in  Eqs.  (96),  (110),  and  (94), 
respectively. 
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In  the  P-E-Y  regime,  Eq.  (130)  becomes  (cos  O'  =“  1) 
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In  the  G-L  regime,  the  analogous  result  is; 
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then  Eq.  (135)  becomes  approximately 
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a 


which  is  formally  identical  to  Eq.  (137).  Thus  Eq.  (139) 
numerical  estimate  of  the  size  of  the  body  to  which  the  G 
proximation  applies;  that  is. 
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provides  a 
L  regime  ap- 
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(141) 


It  also  shows  that  the  G-L  regime  may  be  treated  in  a  manner  similar 
to  the  general  scheme  of  the  P-E-Y  regime.  Therefore,  only  the  P-E-Y 
regime  will  be  examined  in  this  section.  Equation  (135)  may  be  written 
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where  x  is  the  hypersonic  viscous  interaction  parameter.  In  the  pre¬ 
sent  problem,  b/L  is  nondimensional.  In  terms  of  this  representa¬ 
tion,  Eq.  (139)  becomes  Q»  1,  which  would  be  the  condition  for  the 
G-L  regime  approximation.  Values  of  A  and  9^  to  be  used  in  Eq.  (143) 
are  given  in  the  Appendix.  From  Eq.  (92),  one  obtains 


P  V  r 

WWW 

pur 
e  e  e 


n  K(0)  ^ 

TTTT!  1  2 

M 


(145) 


It  can  be  further  shown  that 
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For  very  slender  bodies  of  revolution,  M  -  0(MJ ;  thus,  neglecting 

0(1/M  )  terms, 
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provided  K(0)/I  -  P(l).  Equation  (149)  applies  when  the  effect  of  the 
^w''w*^w^^e“e%  negligible  and  is  of  0(1/M^). 

rne  expression  for  Q  in  Eq.  (143)  contains  e;  thus, 
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These  expressions  can  be  introduced  into  Eq.  (149)  to  yield 
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(152) 


Thus  e  can  be  determined  and  the  hypersonic  viscous  interaction  prob¬ 
lem  is  solved.  To  compute  the  induced  pressure  due  to  viscous  inter¬ 
action,  we  define 
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The  surface  pressure  on  the  body  of  revolution  in  an  inviscid  hyper¬ 
sonic  flow  is  The  induced  pressure  rise  is  due  to  viscous  in¬ 
teraction 
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Therefore 
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VIll.  CONCLUSIONS 

We  have  obtained  a  similar  solution  for  the  case  of  a  hypersonic 
laminar  boundary  layer  on  an  axisynmetric  body.  The  strong  interac¬ 
tion  between  the  boundary  layer  and  the  shock  wave  has  been  considered, 
and  the  analysis  includes  mass  transfer  from  the  surface.  In  the  vis¬ 
cous  flow  regime,  the  boundary  layer  equations  were  transformed  to  a 
set  of  ordinary  differential  equations  employing  appropriate  similarity 
variables.  The  inviscid  region  was  treated  using  the  well-known  simi¬ 
larity  transformation  for  the  small- disturbance  equations.  These  so¬ 
lutions  were  matched  at  the  interface  between  the  two  regions,  requir¬ 
ing  that  the  pressure  and  normal  velocity  be  continuous.  The  matching 
procedure  }  'Ided  an  analytical  expression  for  the  surface  pressure 
and  a  similarity  law  for  the  normal  velocity  at  the  wall. 

We  have  derived  a  set  of  analytic  expressions  that  gives  the  re¬ 
lationship  among  the  induced  pressure  at  the  wall,  normal  injection 
velocity  at  the  wall,  and  other  important  wall  variables,  in  terms  of 
the  parameters  chit  describe  the  strong  interaction  flow. 

A  transformation  appropriate  to  both  the  P-E-Y  regime  and  the  G-L 
regime  has  been  derived  from  Eqs.  (10)  and  (11).  This  transformation 
leads  to  the  general  set  of  equations,  which  can  then  be  applied  to 
the  appropriate  regime. 

Within  the  restrictions  indicated  in  the  present  study,  similar 
solutions  have  been  obtained  for  the  case  of  a  thick  boundary  layer 
on  a  very  slender  body  of  revolution.  However,  when  the  boundary  lay¬ 
er  thickness  and  the  characteristic  body  dimension  are  of  the  same  or- 

(2) 

der  of  magnitude,  it  was  shown,  in  agreement  with  Yashuara's  '  ear¬ 
lier  conclusions,  that  rigorous  similarity  is  possible  only  for  three- 
quarter-power  bodies. 

For  an  axisymmetric  body,  the  similarity  condition  for  the  normal 
velocity  at  the  wall  follows  a  (pvr)  law  rather  than  the  well-known 
(pv)  law  in  the  two-dimensional  case. 

The  restrictive  nature  of  the  similar  solutions  is  adequately  il¬ 
lustrated  in  the  present  study.  In  physical  problems  of  interest,  it 
may  be  necessary  to  deal  with  more  complicated  situations  than  those 
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allowed  b'^  similar  solutions.  We  have  seen  in  Section  VI  that  in  hy¬ 
personic  viscous  interaction  problems,  similar  solution  may  be  appli¬ 
cable  only  in  the  outer  viscous  flow  region  and  may  not  be  applicable 
in  the  inner  viscous  flow  region  (cf.  Eq.  (124d)).  In  such  cases,  we 
strictly  must  obtain  nonsimilar  solutions  of  Eqs.  (38)  to  (40),  under 
the  boundary  conditions  in  Eqs.  (25)  to  (29).  The  scheme  of  calcula¬ 
tions  used  in  Refs.  19  and  20  may  be  adopted  for  this  purpose.  In  the 
present  Memorandum,  the  surface  condition  in  Eq.  (127)  is  adopted  to 
obtain  "approximate"  similar  solutions  of  the  G-L  regime.  Alterna¬ 
tively,  for  the  G-1,  regime,  the  visoous  flow  region  cnay  be  dealt  with 
by  a  composite  layer  approach  in  order  that  all  the  physical  boundary 
conditions  are  rigorously  satisfied.  Stewartson^^^^  has  recently 
treated  the  hypersonic  slender  cone  problems  by  this  composite  viscous 
layer  representation. 
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Appendix 

VALUES  OF  A  AND  0 
■  ■  . . o 

The  similarity  solutions  for  inviscid  hypersonic  flow  over  slen- 
dei  power-law  bodies  have  been  computed  by  Mirels.- In  terms  of 
Mirel's  data  for  axisytnmetric  three-quarter -power  law  bodies,  we  have 


These  values 


Y 

e 

1.4 

.875 

.696 

1.67 

.819 

.634 

of  and 

are  related 

to  0 

o 

and  A  by  the  following 


formulas : 
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2n-l 
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where 
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ny. 


n  -  1  +  nY 


>  0 


3 

4 


43 


REFERENCES 


1.  Li,  T,  Y. ,  and  H.  T.  Nagamatsu,  “Similar  Solutions  of  Compressible 

Boundary  Layer  Equations.**  J.  Aervon,  Sci..  Vol.  22,  No.  9,  Sep- 
•.ember  1955. 

2.  Yasuhara ,  M- ,  "Axi symmetric  Viscous  Flow  Past  Very  Slender  Bodies 

of  Re\ olution. '*  J.  Aeron.  Sci..  Vol.  29,  June  1962,  pp.  667-679. 

3.  Stewartson,  K. ,  '*The  Asymptotic  Boundary  Layer  on  a  Circular  Cyl¬ 

inder  in  Axial  Incompressible  Flow. **  Quart.  Appl.  Math..  Vol.  13, 

No.  '1,  July  1955.  pp.  113-122. 

4.  Glauert,  M.  B.  ,  and  M.  J.  Lighthill,  **The  Axisymmetric  Boundary 

Layer  on  a  Long  Thin  Cylinder,"  Proc.  Roy.  Soc.  London.  Vol.  230, 
1955,  pp.  188-203. 

5.  Mark.  R.  M. ,  '*Laminar  Boundary  Layers  on  Slender  Bodi's  of  Revolu¬ 

tion  in  Axial  Flow,"  GALCIT  Hypersonic.  Wind  Tunnel  L^mo  No.  21, 

July  1954. 

6.  Mangier,  W. ,  "Zusaranenhang  zwischen  ebenen  und  rotationsyymmetrischen 

Greuzschichten  in  Kompressiblen  r lussigkeiten,"  ZAMM.  Vol.  28, 

No,  4,  1948,  pp.  97-103. 

7.  Probstein,  R,  F.,  and  D,  Elliott,  *'The  Transverse  Curvature  Effect 

in  Compresf ■ ble  Axially  Symmetric  Laminar  Boundary  Layer  Effect," 

J.  Aeron.  Sci.,  Vol.  23,  March  1956,  pp.  208-224, 

8.  Steiger,  M.  H. ,  and  M.  H,  Bloom,  **Thick  Boundary  Layers  Over  Slen¬ 

der  Bodies  With  Some  Effects  of  Heat  Transfer,  Mass  Transfer  and 
Pressure  Gradient,"  Int.  J.  Heat  Mass  Transfer.  Vol.  5,  1962, 

9.  Eckert,  E,  R.  G.,  and  P.  J,  Schneider,  "Effect  of  Diffusion  in  an 

Isothermal  Boundary  Layer,**  J.  Aeron.  Sci..  Vol.  23,  Nc.  4, 

April  1956. 

10.  Gazley,  C.,  Jr.,  **Fundamentals  of  Nonequilibrium  Flow  Processes," 

Hypersonic  Gas  Dynamics  Notes,  1958,  pp.  10-13. 

11.  Stewartson,  K. ,  **0n  the  Motion  of  a  Flat  Plate  at  High  Speed  in  a 

Viscous  Compressible  Fluid.  II:  Steady  Motion,"  J.  Aeron.  Sci.. 
Vol.  22,  1955,  p.  303. 

12.  Oguchi,  H. ,  "First  Order  Approach  to  a  Strong  Interaction  Problem 

in  Hj’personic  Flow  Over  an  Insulated  Flat  Plate/*  Aero.  Res. 
Institute,  University  of  Tokyo,  Rep.  330,  1958, 

13.  Li,  T.  Y. ,  and  J.  F.  Gross,  *'On  Transverse  Curvature  Effects  in 
Axisymmetric  Hypersonic  Boundary  Layers,**  AIAA  Journal,  Octobe- 
1964;  The  RAND  Corporation,  P-28/2,  1964. 


•44- 


14.  Li,  T.  V.,  and  J.  F.  Gross,  "Hypersonic  Strong  Viscous  Inter¬ 

action  on  a  Flat  Plate  With  Surface  Mass  Transfer,”  HTFMI  Proc.. 
June  1961<  The  RAND  Corporation,  RM-3000-PR,  March  1963. 

15.  Yasuhara,  M. ,  "On  the  personic  Viscous  Flow  Past  Slender  Bodies 

of  Revolution,  J.  Phys.  Soc..  Japan.  ”ol.  11,  No.  8,  August 
1956  . 

16.  Stewartson,  K.,  "Viscous  Hypersonic  Flow  Past  a  Slender  Cone," 

Phys.  of  Fluids.  Vol.  7,  No.  5,  May  1964,  pp.  667-675. 

17.  Mirels,  H. ,  "Similarity  Solutions  for  Inviscid  Flow  Over  Slender 

Power  Law  and  Related  Bodies,"  ARS  Publication  No.  1111-60,  May 
1960. 

]8.  Li,  T.  Y. ,  and  J.  F.  Gross,  "Hypersonic  Viscous  Interaction  on  a 
Slender  Body  of  Revolution  With  Surface  Mass  Transfer," 
Zeitschrift  fiir  Flugwissenschaf ten.  Vol.  14,  No.  3,  March  1966, 
pp.  129-139^  The  RAND  Corporation,  P-2977-1,  1964. 

19.  Wei,  M.  H. ,  "Asymptotic  Boundary  Layer  Over  a  Slender  Body  of 

Revolution  in  Axial  Compressible  Flow,"  AIAA  Journal.  May  1965, 
Vol.  3,  No.  i,  pp.  809-816. 

20.  Albarete,  L.  M. ,  "The  Effect  of  Fluid  Injection  on  the  Incompres¬ 

sible  Lamiiur  Boundary  Layer  Over  a  Slender  Parabola  of  Revolu¬ 
tion,"  NOLTR  65-226,  February  25,  1966. 


DOCUMENT  CONTROL  DATA 

1 1  ORiGINATING  ACTIVITY 


THE  RAND  CORPORATION 


2a.  REPORT  SECURII Y  CLASSIDCATION 

UNCLASSIFIED _ 

2b.  GROUP 


3.  REPORT  TITLE 


HYPERSONIC  VISCOUS  INTERACTION  ON  A  SLENDER  BODY  OF  REVOLUTION 
WITH  SURFACE  MASS  TRANSFER 


4.  AUTHOR(S)  (Loll  noma,  first  noma, Initial) 


Li,  T.  Y.  and  J.  F.  Gross 


1  5.  REPORT  DATE 


September  1967 

7.  CONTRACT  OR  GRANT  No! 

DAHC15.67.C-0141 
9o  AVAILABILITY/ LIMITATION  NOTICES 


Go.  TOTAL  No.  OF  PAGES 

.S4 

e.  ORIGINATOR'S  REPORT  No. 

RM-4491-ARPA 


6b.  No.  OF  REFS. 

20 


DDC-1 


10.  ABSTRACT 


»b.SPC  >  SORING  AGENCY 


Advanced  Research  Projects  Agency 

II.  KEY  WORDS 


nu  aniilytica'  formulation  of  t  «  problem 
of  hypersonic  viscous  interactiL.;  on  a 
very  slender  body  of  revolution  with  a 
tnick  boundary  layer  and  mass  trails .fer. 
Particular  attention  is  directed  to  tne 
establi  sj.ment  of  tlie  mathematical  restric¬ 
tions  necessary  to  ensure  the  existence 
of  similar  solutions  of  ti.*-  laminar 
boundary-layer  equations  for  this  class 
of  problems.  Yasuhara's  analysis  is  ex¬ 
tended  to  include  surface  mass  transfer 
and  binary  mixture  effects.  It  is  shown 
that  when  the  boundary-layer  thickness 
is  of  the  same  order  of  magnitude  as  the 
characteristic  b~dy  dimension,  rigorous 
similarity  is  possible  only  fo.  three- 
quarter-power-law  bodies.  This  is  in 
agreement  with  Yasuhara's  earlier  con¬ 
clusions.  Similar  solutions  to  the 
boundary-layer  growth  equations  exist 
when  the  transverse  curvature  expression 
is  a  function  of  the  similarity  variabl*. 
The  criterion  for  the  existence  of  simi¬ 
lar  solutions  for  the  problem  of  hyper¬ 
sonic  viscous  interaction  is  that  there 
must  be  a  definite  relationship  between 
body  surface  shape  and  surface  injection 
velocity . 


Ablation 
Mass  transfer 
Hypersonic  flow 
Fluid  mechanics 
Reentry  vehicles 
Aerodynamics 


